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We derive slow-roll conditions for thawing quintessence. We solve the equation of motion of (j> 
for a Taylor expanded potential (up to the quadratic order) in the limit where the equation of state 
w is close to —1 to derive the equation of state as a function of the scale factor. We find that the 
evolution of <f> and hence w are described by only two parameters. The expression for w(a), which 
can be applied to general thawing models, coincides precisely with that derived recently by Dutta 
and Scherrer for hilltop quintessence. The consistency conditions of \w + 1| <C 1 are derived. The 
Oh slow-roll conditions for freezing quintessence are also derived. 
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I. INTRODUCTION 

, The current cosmological observations seems to be consistent with ACDM. The equations of state of dark energy, 
■ w, is close to —1 within 10% or less This implies that even if a scalar field (dubbed "quintessence") plays the 
\ role of dark energy, it should roll down its potential slowly because its kinetic energy density should be much smaller 
I' than its potential. In this situation, as in the case of inflation, it is useful to derive the slow-roll conditions for 
q quintessence. Quintessence models which hardly move in the past and begin to roll down the potential recently are 
^ . called "thawing" models, while "freezing" models move in the opposite ways: they gradually slow down the motion 
0. Our consideration here will be given mostly to thawing models. 

However, as far as we are aware, no such conditions have been derived for general potential. The usual slow-roll 
conditions, (V /V) 2 <C 1 and |V"/V| -C 1, may not be necessary. In fact, in the scalar equation of motion, the 
" acceleration term" , </>, is not necessarily small compared with the " friction term" , 3H </> 0, 0, 01 • The difference 
between inflation and quintessence is that for the former the scalar field always dominates the universe, while not for 
C \ the latter. For slow-roll inflation, the time scale of the scalar field motion is longer than the time scale of the cosmic 
CO . expansion which is determined by the potential. For quintessence, on the other hand, the time scale of the cosmic 
expansion is determined by the matter/radiation and the scalar field. Therefore the potential does not necessarily 
satisfy the usual slow-roll conditions. Note that the situation is not limited to quintessence but is applied to the case 
C^l ' when the scalar fields which are subdominant components in the universe move slowly. Axion, curvaton, and moduli 
, can be such fields. 

If w is close to —1, a functional form of w as a function of the scale factor is necessary to parametrize possible 
deviations from a cosmological constant. The frequently used functional form is the linear approximation of w(a) at 
k*" ■ o=l, the so-called Chevallier-Polarski-Linder parameterization 
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WWar(a) = Wo + W a (l - a). (1) 

Naively, one may expect that this is a good parameterization because it is a Taylor expansion of w(a) at a = 1. 
However, even for a modest redshift z > 2, 1 — a is not so small (1 — a > 0.6) and higher order terms (a— l) 2 , (a— l) 3 . . . 
may no longer be negligible. Recently, it is found that for a particular class of thawing models (hilltop models), the 
equation of motion of the scalar field can be solved explicitly if |1 + w\ -C 1 and w can be written in the closed form 
which is in excellent agreement with numerical solutions Q. It is found that even for dark energy with w ~ —1, its 
equation of state cannot evolve like Eq. ([I]). Unfortunately, however, the analysis is limited to the thawing model 
whose potential has a maximum and appears not to be applied to more general thawing models. We show that this 
is not the case and that the analysis can be applied to more general thawing models with generic initial conditions. 
The functional form of w coincides with that derived in Q. 

In this paper, in Sec. II, we first derive the slow-roll conditions for thawing quintessence models. Then, generalizing 
the analysis in Q, we solve the equation of motion of the scalar field for a Taylor expanded potential (up to the 
quadratic order) in the limit where the equation of state w is close to —1 to derive the equation of state w as a 
function of the scale factor. We also derive the consistency conditions of the approximation. In Sec. Ill, we compare 
the derived w with numerical solutions for several quintessence potentials which unlike Q do not necessarily have the 
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maximum and find fairly good agreement. Sec. IV is devoted to summary. In Appendix ^ the slow-roll conditions 
for freezing quintessence are derived, and in Appendix [B] some calculations which are necessary to derive a equation 
used in the text are given. 

II. SLOW-ROLL THAWING QUINTESSENCE 

Working in units of 8nG = 1, the basic equations in a flat universe are 

4> + 3H<p+V = 0, (2) 

H 2 = (£) =^(PB +/>*), (3) 

H = -\ ((pb+Pb) + (p<i>+P<f>)) = ~\ ((1 + wb)pb + 4> 2 ) , (4) 

where V = dV/dcj), H = a/a is the Hubble parameter with a being the scale factor, pb(pb) is the energy density 
(pressure) of matter/radiation, = <fi 2 /2 + V{(j))(p^ = <f> 2 /2 — V (<fi)) is the scalar field energy density (pressure), and 
Wb is the equation of state of matter/radiation. 



A. Slow-Roll Conditions for Thawing Quintessence 



By slow-roll quintessence we mean a model of quintessence whose kinetic energy density is much smaller than its 
potential, 

^ 2 « V. (5) 

Unlike the case of inflation, we do not require that <j> is smaller than the friction term 3H(j) in Eq. @ since H is not 
determined by the potential alone, but by the matter /radiation along with the scalar field energy density. 

With fixed Wo, slowly rolling thawing models correspond to the equation of state w = p<f>/p4> very close to — 1, so 
that the Hubble friction is not effective and hence cf> is not necessarily small compared with 3H(j> in Eq. (J5J). Slowly 
rolling freezing models correspond to models whose w is not so close to —1 compared with thawing models so that 
the Hubble friction is effective and (f> is smaller than 3Hcj> in Eq. ([2]). 

We derive the slow-roll conditions for thawing quintessence during the matter/radiation dominated epoch. For 
slow-roll conditions for freezing quintessence, see Appendix [A] We first introduce the following function [3| (see also 
1): 

3H<P 

As stated above, for thawing models, (3 is a quantity of 0(1). We assume (3 is an approximately constant in the sense 
|/3| <C H\fi\, and the consistency of the assumption will be checked later. In terms of (3, using Eq. $Z§, <fi is written as 



V 



3{1 + 0)H 

and the slow-roll condition Eq. ([5]) becomes 

V' 2 



(7) 



6H 2 V 



« 1, (8) 



where we have omitted 1 + (3 since it is an 0(1) quantity and introduced the factor of 1/6 so that e coincides with 
the inflationary slow-roll parameter, e = \{V /V) 2 @, if the scalar field dominates the expansion: H 2 ~ V/3. Eq. 
© is a quintessence counterpart of the inflationary slow-roll condition (V /V) 2 -C 1. 

Similar to the case of inflation, the consistency of Eq. ([6]) and Eq. ([2]) should give the second slow-roll condition. 
In fact, from the time derivative of Eq. ([7]) 

' V " V ' (9) 



9(1 + f3) 2 H 2 2(l + /3) 3(l + /3) 2 iT 
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where we have used H/H 2 ~ -3(1 + w B )/2 from Eq. (JSJ) and Eq. g|. On the other hand, from Eq. ((SJ) and Eq. (7}, 
= 3/3-ET0 = -PV'/{1 + /3), and so we obtain 

' 9(l + /3)# 2 2 3(1 +j9)ir - 9(l+/3)# 2 2 1 j 

where we have used $ <g; Hf3. While the left-hand-side of Eq. (JTUJ) is an almost time-independent quantity by 
assumption, the first term in the right-hand-side is a time-dependent quantity in general. Therefore the equality 
holds if the first term is negligible: 



V" 

V ■= 

so that P becomes 



?/:= 3iP ; W<<1 ' 



or the left-hand-side is negligible: 



so that 



0=^, (12) 



« 1, (13) 



V = l(l + WB ). (14) 

The former condition corresponds to the slow-roll thawing models, while the latter corresponds to the slow-roll 
freezing models (see Appendix |A")) . /3 given by Eq. (fl2|) is an approximately constant, which is consistent with 
our assumption. 1 Here the factor 1/3 is introduced in Eq. (JTTJ so that 77 coincides with the inflationary slow-roll 
parameter 77 = V"/V, if H 2 ~ V/3. Eq. (fTl]) is a quintessence counterpart of the inflationary slow-roll condition 
|V"|/V< i- 

Eq. ([8]) and Eq. (jll[) constitute the slow-roll conditions for thawing quintessence during the matter/radiation 
epoch. 2 Moreover once the universe becomes dominated by the scalar field, the two conditions reduce to the usual 
inflationary slow-roll conditions from H 2 ~ V/3. Therefore, these conditions (Eq. (JSj) and Eq. (fTTj) ) are the slow-roll 
conditions for thawing quintessence at all times, both during the matter/radiation era and during the scalar field 
dominated era. Note that since H 2 > V/3, the inflationary slow-roll conditions are sufficient conditions for slow-roll 
thawing quintessence during the matter/radiation era, not necessary conditions. In Fig. 1, the evolution of (3 is shown 
for a thawing quintessence model (V = M 4 (l — cos </>)). The evolution of (3 agrees nicely with Eq. (fTJ 



B. Parametrizing the Equation of State 



Next we derive general solutions of <f> m the limit of |1 + w\ <C 1 and derive w as a function of a. To do so, we first 
note that the Hubble friction term in Eq. ([2]) can be eliminated by the following change of variable Q 

u = (0-^)a 3/2 , (15) 
where fa is an arbitrary constant, which is introduced for later use, and then Eq. becomes 

ti + ^{p B + PlP )u + a 3 / 2 V = 0. (16) 



1 More precisely, $/H(3 = 3(c^ —Wb), where c? s =Pb/pb ls the adiabatic sound speed. Hence the assumption <C H\/3\ is valid if only 
the adiabatic perturbation is considered. 

2 It is easily found that even if multiple scalar fields each with the same potential are introduced, the assisted dynamics lid does not 
occur: e and r\ remain the same. 
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FIG. 1: (3 as a function of a for thawing quintessence model with the axion-like potential V = M (1 — cos (f>). The dotted lines 
are j3 = 2/3, 1/2, respectively. 

We assume a universe consisting of matter and quintessence with w ~ —1. Then the pressure is well approximated 
by a constant: pg + p<p ~ p^ ~ —p<f,oi where p^o is the nearly constant density contributed by the quintessence in the 
limit w ~ —1. Eq. (|16p then becomes 

u - ^ P<f>0 u + a 3 / 2 V = 0. (17) 

Since wc consider a slow- roll scalar field, the potential may be generally expanded around some value (f>i, which we 
identify with the initial value, in the form (up to the quadratic order) 

V{4>) = Vfa) + V'ttM - + \v"{4>i)(<j> - 4) 2 - (18) 

We will check the consistency of the expansion later. Substituting the expansion Eq. (fTS)) into Eq. (fT7|) and taking 
P<t>o = V{4>i) gives 

u + (v"(<k) ~ \v{<t>i)\ u - -F'(^)« 3/2 . (19) 

Here the source term in the right hand side of Eq. (fTT))) appears since we consider the general Taylor expansion of V 
around the initial value <pi in contrast with Q where the Taylor expansion of V around its maximum is considered 
and hence V term in Eq. (fT5)l is absent. 
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Being consistent with \w + 1| -C 1, we assume a(t) is well approximated by its value in the ACDM model which is 
given by 



o(t) 



1 - Q 



00 



SI 



1/3 



sinh 2 / 3 (i/i A ), 



(20) 



where SIaq is the present-day value of density parameter of quintessence, a — 1 at present, and £a is dchncd as 



tA 



(21) 



Introducing 



k = V(3/4)V(&)-V™(&), (22) 
the general solution of Eq. (fTT)|) is obtained by the use of Green function method in the form (if kt\ ^ 1) 

'i-n 0O v'{<t>i)tl 



i(t) = Asmh(kt) + Bcosh(fci) + 



S1 0O kH\ - 1 



sinh(i/i A ), 



(23) 



where A and B are constants. kt A — 1 corresponds to V"(4>i) = 0, which will be treated separately later. As an 
initial condition, we take that <f> = <f>i and <j) = 4>i at t = ti. Then, we obtain 



<Kt) - <t>i 



sinh (^ /t a) 
kt\ sinh(t/<A) 



sinh(fct) cosh(fcti) 



V'{<l>i) 



V"(4>i) 



(coth(fi/<A) — kt A tanh(fcf,)) + t A 4>i 



cosh(fci) sinh(fctj) <^ (coth(i;/t A ) - kt A coth(fctj)) + t A <& 

1 V"(0i) 



(24) 



However, as shown in Appendix [Bl as long as <Zj <C 1 (or ti <C £o), the solution can be well approximated by that with 
U=0 



0(*) = & 



sinh(fct) 



V"{<t>i) \kt A sinh(t/t A ) 
Taking ~ p^o — V'(^i), using Eq. (|25p the equation of state is given by 



l + w{t) 



<j) 2 3/ V(<&) \ /W A cosh(fci)sinh(i/i A ) - sinh(fci) cosh(i/i A ) 



Vfa) 4 V*tAV"(0i) 



sinh (t/tA) 



(25) 



(26) 



As done in Q, we normalize the expression to the present-day value of w, Wq, and rewrite masa function of the scale 
factor using Eq. (|20j) . Normalize Eq. (f2"6"| to the present-day value, 



1 + w(a) = (1 + w )a 



-3 



if cosh(fei(a)) - F(o) sinh(fct(a)) 
if cosh(fcio) — ^0o^ 2 sinh(fcio) 



(27) 



where if = kt A and -F(a) is the inverse square root of the fractional energy density corresponding to a cosmological 
constant and they are given by 



/ 4 V" 
R = kt a = 4 / 1 - - 



3 Vfa) ' 



F(a) = yjl + (Sl^ 1 - l) a -3. 



t(a) can be derived from Eq. ([20]) so that 



kt(a) = K sinh 1 J — 



a 3 SI 



</)f) 



S> 



if In 



i 3 SI 



00 



00 



1 - SI 



</<() 



(1 + F(a)) 



(28) 
(29) 

(30) 
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FIG. 2: G(a = 1, K) as a function of K for ft^ = 0.74. 



Then Eq. (|2"Tj) can be written as 



, * ri , x 3(k-d / (g - f(«))(^(q) + i)* + + WXW - i) g 

1 + w(a) = (1 + woja v ' 775 TF> TB TF> 



(31) 



Remarkably, the expression Eq. (|3 1 [) formally coincides with that of 0] (Eq. (31) in 0]) where the expression is 
derived for hilltop quintessence. However, the definition of K is different: K = y/l — (4/3)V"/V in Q where V and 
V" are evaluated at the maximum of V, while our (Eq. (|28[) ) is evaluated at the initial value Moreover, we 
derive this expression for a general Taylor expanded potential around <f)i- Hence Eq. (|3ip is not limited to hilltop 
quintessence but can be applied to a much wider class of slow-roll quintessence, ft should be noted, however, that 
since w(a) is an increasing function of a, which is evident from the derivation, Eq. (|31|) can be applied only to thawing 
model: model with growing w. 

We also note that we can consider the case of K 2 < (or V" > (3/4)V) by the replacement: K — iK' where 
K' = y/{4/3){V"/V) - 1. Using Eq. (001), w(a) corresponding to Eq. {TQ is give by 

( K' cos (k' sinh- 1 y/jlzjfc) - F{a) sin (k 1 sinh" 



1 + w(a) = (1 + w )a 



-3 



^ K> cos {K> sinh- 1 y 7 ^) - n~l /2 sin (l<> sMT 




(32) 
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C. Consistency Conditions 



Now, we check the consistency of the expansion Eq. (|18p and the assumption |1 + w\ <C 1. The expansion Eq. (Tig]) 
is consistent if V(<j>i) > \V - fc)] > | (1/2) V" (&)(</> - <fo) 2 |. Using Eq. this is satisfied if 



r := 
G(a,K) :- 



sinh(fci) 



A"sinh(i/i A ) 



< 1, 



(33) 
(34) 



where we have introduced T first introduced (without the absolute value) in |Tl| . Note that G(a,K) is an increasing 
function of a. Using Eq. (j3"0")l . G(a = 1, K) is given by 



G(a = 1,K) 



1 



2K v i - n 



16 



(K-l)/2 



((^o /2 + 1)* - (^0 



1/2 _ ^if 



- 1 



(35) 



In Fig. 2, G(a = 1,K) is plotted for f^o = 0.74. It is smaller than unity and the inequality Eq. (|34[) is satisfied if 
\K — 1| < 1. In terms of the potential, this implies 



V{<t>i) 



< 1. 



(36) 
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Since the epoch of the initial conditions is arbitrary as long as a* <C 1, 4>i hi Eq. (j33|) . Eq. 
with <p. Then we finally obtain the consistency conditions of slow-roll thawing quintessence: 



r = 



v(4>)v"(<j>) 



v>(cf>y 



= o(i), 



V"{<t>) 



< 1. 



36} may be replaced 

(37) 
(38) 



Note that from (V /V) 2 = T-^V" /V\, these two conditions imply (V /V) 2 < 1. Since 3H 2 > 
ter/radiation epoch, \V" /V\ < 1 and (V 1 /V) 2 < 1 are consistent with the slow-roll conditions, Eq. {£ 
which ensures |1 + w\ <C 1. 



V during mat- 
} and Eq. {IT}, 



D. K = 1 Case 



Lastly, for completeness, we consider the case of K = 1 (or V"(0j) = 0) Q. In this case, the general solution of 
Eq. {19} with <j) = (f) t at t = is 



0(t) = ^ + 



2V'{ 



3V(&) 



The equation of state is 



1 + w 



3 V Vfa) 

( 

(l + wo) 



/ sinh(fct) cosh(fci) — fci x 2 
sinh 2 (fci) 



F(a) 



hi 



V 



Therefore, the Taylor expansion is consistent if 



O 



-1/2 l-Qj 



hi 



i+ni 



< 1, 



(39) 



(40) 



(41) 



which is again consistent with the slow-roll condition, Eq. {8}, and hence ensures |1 + w\ -C 1. Eq. (140(1 is identical 

to the expression derived for a linear potential in Q . It is easily found that Eq. (|3"Tj) reduces to Eq. (14T)} in the limit 
K — > 1, confirming the result by Q. 

In Fig. 3, ui(a) given by Eq. (|3"Tj) or by Eq. {40} is shown for several K . 



III. COMPARISON 



A. Comparison with Numerical Solutions 



We compare the slow- roll prediction of w(a) (Eq. {3T}) with numerical solutions for several models and evaluate 
the accuracy of Eq. (f3"T} . We consider the following three examples: 

(a) the pseudo Nambu-Goldstone boson (axion-like) model [12 ]: 

V(0) = Af 4 (1 - cos(0//)) , (42) 

(b ) logarithmic potential [l3[ : 



y(^) = M 4 iog(0), 



(43) 
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FIG. 4: w(a) for (a) the axion-like potential, V — M 4 [l — cos0], for (b) the logarithmic potential, V = M 4 log0, and for 
(c) the quadratic potential, V = m 2 4> 2 /2. The solid (black) curve is the numerical solution, the dotted (blue) curve gives our 
approximation, and the dashed (red) curve gives the approximation in Q. 



(c) quadratic potential [Til ]: 

V(4>) = im 2 2 , (44) 

where M, f and m are constants. The first example is considered in 0], and the second example corresponds to the 
potential without maximum/minimum, and the third example corresponds to the concave potential V" > so that 
K < 1. We fix fi^o = 0.74 and take / = 1 in the reduced Planck units and choose <fii so that Wq ~ —0.9. The results 
are shown in Fig. 2J 

For all cases, we find fairly good agreement with the numerical solutions: For case (a), the relative error (the 
difference between our approximation and the numerical solution), \5w/w\, is less than 0.3% while it is less than 0.1% 
for the approximation by For case (b), \5w/w\ < 0.3% and for case (c) it is less than 0.7%. Note that the potential 
does not have a maximum for the latter two cases and the approximation of Q is no longer available. To check the 
slow-roll conditions Eq. §7} and Eq. ([38]), we compute T, V" /V at T = 1.33, 0.77, 0.50; \V"/V\ = 0.31, 0.28, 0.22, 
respectively. 3 



B. Comparison with Other Parametrizations 



Finally we compare our parametrization with other parametrizations of w(a). 

The most frequently used functional form is the linear approximation of w(a) at a = 1, the so-called Chevallier- 
Polarski-Linder parameterization, wn ncar (a), Eq. |T]) Q. 

Another parametrization of w(a) closely related our approach is that by Crittenden et al.@. Instead of expanding 
the potential, they expanded the slow-roll parameter around cf>o in linear order: 

V 

K W = q - «y = + «i(</> - <hh ( 45 ) 
Resulting w, denoted as w cmp (a), is written as Q 



l + w cmp (a) = i«^; 2Kl/3 a- 2 -F( a )-( 4 -+ 6 )/ 3 



3 

(1 + w )%^ 1+3)/3 a- 2 ^F(a)-^ +e ^ 3 , (46) 



3 It is to be noted that the slow-roll conditions Eq. 11371 and Eq. (13S I t arc required for the consistency of the solution Eq. (1 3 IP and define 
the range of validity of the solution; otherwise the expansion of the potential Eq. I I18H and |1 + w\ <C 1 cannot be trusted. However, this 
does not imply that Eq. H31H cannot be used for \V" /V\ S> 1 (or K 3> 1), but rather simply that we cannot trust such an extrapolation. 
Interestingly, for axion case (and other hilltop quintessence model), it is shown that approximation Eq. I|31fl with K evaluated at the 
maximum of the potential agrees excellently with the numerical solution even for K = 4 0] . 
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FIG. 5: Top: Comparison between our parametrization w(a) (Eq. Q31| l) (solid curve) and a linear parametrization Wii near (a) 
(dotted curve) and «; cmp (a) 0| (dashed curve) for K = 0.8, 1, 1.2 from left to right. Bottom: the relative error between wimcar(ct) 
(dotted) or w cmp (a) (dashed) and our parametrization: Sw/w — (iuii ncar / cm p — w(a))/w(a). K — 0.8, 1, 1.2 from bottom to top. 

where in the last line we have normalized the equation of state to the present-day value. 

In the upper panel of Fig. [5j three parametrizations of the equation of state normalized to the present-day value 
(and the first derivative) are shown. In the lower panel, the relative error is shown. It can be seen that the difference 
between the linear parametrization wn nea , T (a) and ours is less 2% for a > 0.5 (or z < 1), but it can be as large as 
6% for smaller a. On the other hand, the difference between w cin p(a) and ours is at most less than 1%. Hence, as 
far as the goodness of fit is concerned, there is no difference between them. However, while for w cmp (a) Ki is related 
(roughly) to the first and the second derivative of the potential, for our parametrization Eq. (pTTj) K is directly related 
to the curvature of the potential. 

It should be stressed that our w(a) Eq. pTj) is not a fitting function particularly designed to match the numerical 
solutions, but rather a function derived by solving the equation of motion of 4>- Our results demonstrate that only 
a slight change of the definition of K in [3] as in Eq. (f!?5|) greatly extends the applicability of Eq. (pH) . Eq. (pIT]) 
can be used not only for hilltop quintessence but also for other quintessence with a more general potential without 
maximum. Conversely, we propose that the parametrization of the equation of state of the form Eq. (|3ip with two free 
parameters (wq, K) may be used to fit the cosmological data. It fits better than the commonly used linear equation 
of state, Wiinear(o) = Wo + W a (l — a), and more importantly the meaning of the parameter K is clear: the curvature 
of the potential. 
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IV. SUMMARY 

We have derived slow-roll conditions for thawing quintessence models, Eq. ([8]) and Eq. (fTTj) . We have also solved 
the equation of motion of the slow-roll thawing quintessence and obtained the equation of state as a function of the 
scale factor w(a), Eq. (f3"Tj) . which involves only two parameters. We have derived the consistency conditions of the 
approximation, Eq. (|3T[) and Eq. (|38|) . which are consistent with the slow-roll conditions. We have found that only 
a slight change of the definition of K in 0] greatly extends the applicability of their w(a). We have shown that our 
w(a) agrees fairly well with the numerical solutions for several thawing models and found that our w(a) is in general 
not fit by a linear evolution in a as emphasized by Q ■ 

It would be desirable to have useful approximation of w(a) for freezing quintessence models and to obtain the 
unified expression for w(a). However, to do so, the different approach is needed, since the equation of state can 
be significantly different from —1 during matter/radiation epoch. We have derived slow-roll conditions for freezing 
quintessence models, Eq. and Eq. (|A4|) . 

It would also be interesting to extend the slow-roll conditions to quintessence with non-minimal coupling with 
gravity (extended quintessence) [l5| by extending the results for non-minimally coupled inflaton(s]_ 1(| , which could 
provide conditions for "tracking without tracking" to solve the coincidence problem dynamically [17 1. 
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APPENDIX A: SLOW-ROLL CONDITIONS FOR FREEZING QUINTESSENCE 

We derive the slow-roll conditions for freezing quintessence during the matter/radiation epoch. For slowly rolling 
freezing models, although the kinetic energy density is smaller than the potential Eq. (0), w is not so close to — 1. 
Then, compared with slowly rolling thawing models, the Hubble friction is effective and so in Eq. ([2]) 

« 1, (Al) 



3H<f> 

and the usual slow-roll equation of motion is obtained: 

3H <j> + V' = 0. (A2) 

Using Eq. (|A2[) . the slow- roll condition Eq. ([5]) becomes again e <C 1, (Eq. (JSJ). The time derivative of Eq. (|A2[) 

gives 

Therefore, the assumption Eq. (|A1[) is consistent if 

V" 3 

V = 3IP = 2 (1+Wb) ' (A4) 

which coincides with Eq. (fl"4")) . Eq. |8]) and Eq. (|A4[) constitute the slow-roll conditions for freezing models. In Fig. 
[6l the evolution of w, (3 and r\ for a freezing model (V = M 4 </> -1 / 4 ) is shown. Although rj deviates from Eq. (|A4|) to 
the extent that (3 deviates from zero (or w deviates from —1), we find reasonable agreement with Eq. (|A4j) . 

APPENDIX B: REDUCTION OF EQ. fl24]) TO EQ. (p5j) 

The general solution of Eq. (TT5)) with (j) — 4>i and = at t = t t is given by Eq. (|2"4"|) . Using Eq. (|2"D|) and Eq. 
(pOl) . we obtain 

e kt > = ( T ^2_) A/2 af/ 2 (F(a l ) + l)^. (Bl) 
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FIG. 6: w, j3 and 77 as a function of a for a freezing quintessence model with the potential V = M 4 -1 ^ 4 . The dotted lines are 
77 = 2, 3/2, respectively. 



For dj <C 1, since F(fli) ~ \/^m ~ fr° m Eq. (f2U|) . we obtain 



sinh(fctj) ~ K> 



1 - n 



<;»() 



cosh(fc^) ~ l + 0(af). 



We also note that 



sinh(i !: /< A ) = l/ r^TT » cosh(^/i A ) = F(aj) sinh(ij/*A) 
y 1 — S2^o 



(B2) 



(B3) 



Hence, in the limit of a, -C 1, only the term proportional to smh(kt) survives in Eq. {H} (up to C(a- /2 )). Thus, Eq. 
(|2il) is reduced to Eq. ([25]) . 
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